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Abstract — In this paper, we mainly focus on the 
transmission of dengue fever with the aim of analyzing and 
comparing logistic and exponential curve that occurred in 
Lahore (Pakistan), in the year 2010—2011. We obtain the 
force of infection, A, from the real data of Lahore, and then 
computed the basic reproduction number, R . Based on the 
analysis of the behaviour of R , we show that the disease is 
epidemic in the region. We also compare the values of 
logistic fit and exponential approximation to find their 
relative percent error. 
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I. Introduction 

Dengue fever, a mosquito-borne viral infection found in 
tropical and subtropical regions worldwide. The World Health 
Organization estimates that today more than 2.5 billion people 
are at high risk for Dengue fever [18]. Currently the disease is 
endemic in more than 100 countries, and is spreading 
throughout the world. Dengue disease is caused by four 
antigenic-ally related but distinct serotypes namely DEN-1, 
DEN-2, DEN- 3 and DEN-4 in which only DEN-2 and DEN-3 
are mostly found in tropical countries [5]. It is transmitted to 
the human body by biting an infected female Aedes Aegypti or 
Aedes Albopictus mosquitoes as the primary mosquito vector. 
Infection by any single type apparently produces permanent 
immunity to that particular serotype to which the patient is 
exposed, but only temporary cross immunity to the other three 
serotypes (Cunha, 2007; WHO, 2002). Infection with one of 
these viruses i.e., DEN-1, DEN-2, DEN-3 or DEN-4, 
characteristically results in fever, headache and rash. The 
clinical spectrum can vary, however, from asymptomatic to 
more severe infections with bleeding and shock. A second 
infection by this virus can result in a more virulent form of the 
disease known as dengue hemorrhagic fever (DHF). 

In twentieth century the disease was common in Asia and 
Pacific [17]. In Asia, the first confirmed outbreak of DHF 



occurred in Philippines and Thailand in the year 1950. 
However, in the mid of the year 1970, the disease spread 
throughout South East Asia. Thereafter in the next 20 years, 
dengue transmission is epidemic every 3-5 years in hyper 
endemic areas [15]. 

In Pakistan, the first confirmed dengue haemorrhagic 
fever outbreak was due to serotype DEN-2 reported in June 
1994 by Aga Khan University Hospital (AKUH) [4]. 
Thereafter, a large number of DHF cases have been reported 
from different parts of the country. Akram et al. reported that 
two serotypes namely DEN-1 and DEN-2 were founded in the 
sera of children with undifferentiated fever [1]. In 1998, an 
outbreak of dengue fever in the Baluchistan province was 
occurred due to co-circulation of DEN-1 and DEN-2 [14]. In 
year 2005, a number of DHF cases due to serotype DEN-3 
were reported in different hospitals of Karachi [8]. Thereafter, 
in year 2006, about 3,640 cases of DHF due to serotypes 
(DEN-2 and DEN-3) were admitted to several hospitals in the 
country [11]. 

The dengue disease transmission cycle starts with a dengue 
infected host body. The virus circulating in the blood of these 
infected bodies in the viremic period (4 to 7 days) called 
intrinsic incubation period [6, 8]. During this viremic period, an 
uninfected female Aedes aegypti mosquito bites a human body 
and ingests blood that contains dengue virus. Then an extrinsic 
incubation period begins within the mosquito (vector). 
Extrinsic incubation period is about 8 to 12 days. After that the 
salivary glands of the mosquito become infected and the virus 
is transmitted when the infected mosquito bites and injects the 
salivary fluid into the wound of the human body. 

The life span of Aedes aegypti mosquito is about fourteen 
days. Infected female mosquitoes never recover from the 
infection since their infective period ends with their death [6]. 
As no vaccine presently exists, the only method of controlling 
or preventing dengue and DHF is to combat or eradicate the 
mosquito vectors [18]. Relative humidity and temperature of 
the region affects the development of the mosquitoes as well as 
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dengue viral development. As we see in Pakistan, dengue cases 
are generally at peak in June, July and August. To control the 
dengue and dengue hemorrhagic fever effectively, we should 
understand the dynamics of the disease transmission and take 
into account all of the relevant details. 

In this paper, we used a dengue fever transmission model 
that was earlier proposed by Pinho et al. in 2010. The human 
population is divided into susceptible, exposed, infected and 
recovered classes, while the vector population is divided into 
aquatic, susceptible, exposed and infected classes. The 
proposed model is tested on real epidemic data of Lahore 
(Pakistan). The model shows that there is a reproductive ratio 
Ro that conceptualizes the rate of spread of a dengue disease 
and determines a threshold: whenever Ro < 1, a typical 
infective give rise, on average, to less than one secondary 
infection, and the disease will die out. Otherwise when Ro > 1, 
the disease will persist in the population [10]. 

II. A MATHEMATICAL MODEL FOR DENGUE 
FEVER IN A VIRGIN ENVIRONMENT 

In this paper, we use a Dengue fever model proposed by 
Pinho et al. in 2010 [16]. In this model, host population N H (t) 
is divided into four compartments, one is known as susceptible 
(a class of human which is exposed to be infected), denoted by 
S H (t), second one is exposed class which is denoted by E H (t), 
third is infected class which is denoted by I H (t), and forth a 
recovered class (those individuals which are from diseases or 
removed), denoted by R H (t). Vector (mosquitoes) population 
N v (t) is divided into four compartments, one is called aquatic 
class, denoted by A v (t), second one is a class that potentially 
infected by dengue virus S v (t), third class is exposed class 
denoted by E v (t) and mosquitoes that were infected with 
dengue virus I v (t). 

The Pinho model is a system of eight ordinary 
differential equations. The block diagram of the model is given 
in Figure (1). Using Block diagram in Figure (1), the dynamical 
system for host (human) and vector (mosquitoes) population 
are described by the following two systems of ordinary 
differential equations: 



(1) 



dS H (t) 

dt 
dE H (t) 
d(t) 

dl„(t) 

dt 
dR H (t) 

dt 

where X is the birth and death rates of human population, p. 
is the mortality rate of vector population, a is the mean number 
of bites per day per vector per host, [3m is the effective contact 
rate, ca and cm are the control effort rates, C is the vector 
carrying capacity, y is the inverse of the duration of host 
viraemia, \x is the mortality rate of the vectors, (j.a is the average 
aquatic mortality rate, ym is the mean aquatic transition rate, k 
is the hatched from all eggs, 8 is the mean oviposition rate, 8e 
(extrinsic) and 9i (intrinsic) incubation rates of vector and host 
populations respectively. 
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Figure 1 . Block Diagram of SEIR Dengue Fever Model 

=kS(l-A v (t)/C)N v (t)-(r m +M a +c a )Mt) 

=r m Mt)-^ m s v (t)i H (t)/N H (t)-( M +c m )s v (t) 

= aJ3 m S v (t)I H (t)/N H (t)-(fi+e e +c m )E v (t) 



(2) 



= e e Ey(t)-(M + C m )Iy(t) 



d\(t) 
dt 

dt 
dEyjt) 
dt 

dlyjt) 

dt 

with two constraints 

N H (t) = S H (0 + E H {t) + I H (t) + R H {t) 

and 

N v (?) = \(t) + S v (0 + E v {t) + I v {t) 

III. LOGISTIC CURVE FITTING 

Since, the Lahore dengue epidemic data is of logistic type, 
therefore using curve fitting method, the logistic curve can be 
described by the equation: 



I H (t) = 



L 



l + expC/^+ZV) 



(3) 



Where I H (t) is the total number of infected human at time t, 
L is the carrying capacity (upper limit) of the population, and 
(3 and (3[ are variables that affect the shape and orientation of 
the logistic curve [13]. 

Using Lahore dengue epidemic data points from table (1) 
we obtained values of (3 = 1.911 and Pi = -0.1348. The value 
of carrying capacity was chosen as L = 19,220, using a 
combination of standard fitting techniques and visual 
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inspection [3]. The resulting logistic equation, fitted for the 
Lahore epidemic data, is 



19220 



l + exp(l .91 1-0. 13480 



(4) 



TABLE I. 



CUMULATIVE NUMBER OF INFECTED CASES FOR 
DIFFERENT HOSPITAL OF LAHORE 



Days 


Cumulative number of 
infected cases 





2476 


12 


3577 


24 


9585 


40 


11966 


50 


14981 


60 


17605 


71 


19211 



IV. THE FORCE OF INFECTION AND EXPONENTIAL 
CURVE 

"A" (the force of infection) is defined as a measure of the 
risk of a susceptible person to become infected, per unit of time 
t [12]. For small values of t, initially in the beginning of an 
epidemic, the infected number of accumulative cases can be 
vary as 



I H (t)ccexp(At) 



(5) 



V. 



LOGISTIC CURVE and EXPONENTIAL 
APPROXIMATION 



The total number of infected human at a given time t is 
estimated accurately using a logistic fit and it is also 
approximately equal to exponential curve as shown [13] 



l + exp(/? +/?,0 



2 
1.8 
1.6 
1.4 
1.2 

1 

0.8 
0.6 
0.4 
0.2 



oc exp(-/?,0, for sma11 1 



(6) 



10 20 30 40 50 60 



Figure 2. Cumulative number of infectious verses time (in days). 

www.ijtel.org 



TABLE II. RELATIVE PERCENT ERROR BETWEEN LOGISTIC FIT 
DATA AND EXPONENTIAL APPROXIMATION 



Days 


Actual logistic fit 


Approximation 


Relative % error 


1 


T7QO C 


Zoj j.j 


1 .ozo 


Z 


J 1 lo.y 


JZ4Z.Z 


j.yjj 


i 
J 


1AQH 1 
J4o /.J 


3710.0 


o.ioo 


4 


3888.8 


4245.4 


9.169 




4323 7 


4858 1 


1Z.JJ7 


6 


4792.1 


5559.1 


16.005 


7 


5293.2 


6361.4 


20.180 


8 


5825.5 


7279.4 


24.957 


9 


6386.8 


8329.8 


30.422 


10 


6974.0 


9531.9 


36.678 




Figure 3 . Data 1 shows logistic curve and data 2 shows exponentional curve 

Figure (3) and Table (2) show the visual and statistical 
comparisons of the logistic and exponential approximations. 
From equation (5) and equation (6) it can be seen that A = - fa = 
0.1348 [3]. 




Jan FebMarchApril May June July Aug Sep Oct Nov Dec 

Figure 4. Monthly temperature of Lahore (Pakistan) 
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We collected metrological data from the metrological 
department of Punjab province during the epidemic year 2010- 
11. Data include temperature (C°) of Lahore. The average 
yearly temperature of Lahore is 24.5 C . Graph of temperature 
of Lahore against months is shown in figure (4). 

We also collected dengue fever data during the year 2010- 
11 from directorate of general health Punjab, which is used in 
the estimation of basic reproductive number R . 

VI. ESTIMATION OF BASIC REPRODUCTION 
NUMBER R 

The basic reproduction number R is defined as "the 
number of secondary cases when an infected person is 
introduced in a totally human susceptible population". 
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(7) 



The parameters which are used in the estimation of R are 
taken from Pinho parameter table (3). The force of infection 
"A" is computed from equation (5) and equation (6). "A" 
influence the basic reproduction number R [2]. 



TABLE III. 



PINHO PARAMETER VALUES 



Parameter 


Biological name 


Parameter 
ranges 


Mean 
value 


e e 


Extrinsic incubation rate 


0.02-0.2 


0.11 


St 


Intrinsic incubation rate 


0.083-0.17 


0.126 




Vector mortality rate 


0.02-0.09 


0.055 


A 


Human mortality rate 


0.0143-0.0167 


0.015 


r 


Human recovery rate 


0.083-0.25 


0.166 


Cm 


Control effort rates 


0-1 


0.5 


s 


Average oviposition rate 


0-11.2 


5.6 


Ym 


Average aquatic transition rate 


0-0.19 


0.095 


k 


Fraction of female mosquito 
hatched from all eggs 


0-1 


0.5 


c a or c m 


Control effort rates 


0-1 


0.5 


a 


Average bite per mosquito per 
day 


0-1 


0.5 


pm 


Effective contact rates 


0.75 


0.75 



In Table (4), we estimated R by using Pinho et al. 2010 



[16]. 



TABLE IV. ESTIMATION OF R FOR DENGUE EPIDEMIC (EP) 
DATA OF LAHORE HOSPITALS (LH) (PAKISTAN) M (PINHO ET 
AL.(2010)) 



EP 


A 


X 


r 




Cm 


e. 


e e 




LH 
data 


0.13 


0.02 


0.17 


0.06 


0.5 


0.13 


0.11 


2.2 



VII. DISCUSSION and CONCLUSION 

Dengue has induced havoc impact on tropical and 
subtropical regions of the world. For the last few years it has 
badly affected the some areas of our country Pakistan 



especially Punjab province. This area remained under sever 
attack of dengue fever in 2010— 201 1. 

In this research study, we monitored the scheme of spread 
of disease in Punjab province by using the proposed 
mathematical model of Pinho et al. (2010). In the beginning of 
epidemic due to correlated symptoms of dengue with flu, very 
few cases were registered but later on due to efficient effort of 
Punjab government lot of data were there in hospitals for our 
research study. We used some key parameters (Pinho 
parameters) which presents the spread of disease, other 
parameters are directly calculated from Punjab data e.g., Po = 
1.91 1, Pi = -0.1348 and L = 19220. Using these parameters and 
force of infection (A = -p^. With these values we calculated the 
basic reproduction number R . Value of R describes the spread 
of dengue disease in Punjab. By using technique of this 
mathematical model one can predict the peak number of 
infection per day and the time at which it happened. 

From our calculation of R , we see that dengue is epidemic 
in Punjab (R >1) during monsoon period. Model suggests that 
spread of disease can be controlled by reducing the oviposition 
rat (§) and exposed humans E H (t). This can be done in very 
low expanses like reducing the places of stagnant water, 
reducing the migration of infected people I H (t), by using 
mosquito's nets. 
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